Polyakov Loop and Gluon Quasiparticles in Yang-Mills Thermodynamics 
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We study the interpretation of Lattice data about the thermodynamics of the deconfinement 
phase of 5(7(3) Yang-MiUs theory, in terms of gluon quasiparticles propagating in a background 
of a Polyakov loop. A potential for the Polyakov loop, inspired by the strong coupling expansion 
of the QCD action, is introduced; the Polyakov loop is coupled to tranverse gluon quasiparticles 
by means of a gas-like effective potential. This study is useful to identify the effective degrees of 
freedom propagating in the gluon medium above the critical temperature. A main general finding 
is that a dominant part of the phase transition dynamics is accounted for by the Polyakov loop dy- 
namics, hence the thermodynamics can be described without the need for diverging or exponentially 
increasing quasiparticle masses as T — > Tc, at variance respect to standard quasiparticle models. 

PACS numbers: 12.38.Aw,12.38.Mh 
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I. INTRODUCTION 

Since the seminal paper by Yang and Mills the 
study of nonabelian gauge theories has attracted a lot 
of interest in several sectors of physics. Focusing on the 
gauge theory of strong interactions (QCD), thanks to im- 
provement of computer facilities which allovif to perform 
huge simulations of Yang-Mills theories, as well as to the 
possibility to create in laboratories extremely hot envi- 
ronments by means of heavy ion collisions, the interest in 
understanding the thermodynamical properties of strong 
nonabelian mediums has noticeably increased. 

Lattice simulations of SU{3) Yang-Mills theories have 
been performed, see for example Refs. 0411 ■ All of these 
studies agree about the onset of the deconfinenemt phase 
transition at T = w 270 MeV. Below T^, the thermo- 
dynamics of the Yang-Mills theory is dominated by the 
lowest lying glueballs. For what concerns the high tem- 
perature phase, the picture is not so clear. As a matter 
of fact, one of the most intriguing aspects is the nonper- 
turbative nature of the gluon medium above the critical 
temperature 0. For example, in Q it is pointed out 
that the perturbative regime is realized for temperatures 
T > lOOTc (more concretely, at T = lOOTc both pres- 
sure and energy density are below the Stefan-Boltzmann 
limit, which corresponds to the ideal gas case, by about 
10%). This makes the identification of the correct de- 
grees of freedom of the gluon plasma, in proximity as well 
as well beyond the critical temperature, a very compli- 
cated task. Resummation schemes have been proposed, 
based for example on the Hard Thermal Loop (HTL) ap- 
proach [rl-lioj. More recently within a NNLO in the HTL 
perturbation theory ll| it has been shown that an agree- 
ment with the lattice results can be obtained also for the 
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trace anomaly down to T ~ 2 Tc but with large uncer- 
tainties. At these high temperatures, the HTL approach 
motivates and justifies a picture of weakly interacting 
quasi-particles, as determined by the HTL propagators. 

Within this theoretical framework, the quasiparticle 
approach to the thermodynamics of QCD has attracted 
a discreet interest recently, see for example Refs. [l^ - [27j 
and references therein. In such an approach, one iden- 
tifies the degrees of freedom of the deconfinement phase 
with transverse gluons; the strong interaction in such a 
non perturbative regime is taken into account through 
an effective temperature-dependent mass for the gluons 
themselves. This theoretical approach is similar to the 
one introduced by Landau \2§\ to describe strongly cor- 
related systems. Generally speaking, it is not possible to 
compute the self-energy of gluons exactly in the range of 
temperature of interest; as an obvious consequence, the 
location of the poles of the propagator in the complex 
momentum plane, hence (roughly speaking) the gluon 
mass, is an unknown function. For this reason, one usu- 
ally assumes an analytic dependence of the gluon mass 
on the temperature, leaving few free parameters which 
are then fixed by fitting the thermodynamical data of 
Lattice simulations. 

The advantage of such an approach is, at least, 
twofolds. Firstly, it is of a theoretical interest in itself 
to understand which are the effective degrees of freedom 
of the Yang-Mills theory at finite temperature and in par- 
ticular the evolution of non perturbative effects with tem- 
perature. Furthermore, once a microscopic description of 
the plasma is identified, it is also possible to include the 
implied dynamics into a transport theory capable to di- 
rectly simulate the expanding fireball produced in heavy 
ion collisions computing the collective properties, as well 
as the chemical composition of the fireball as a function 
of time [H]. 

In this article we extend the quasiparticle picture of 
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the finite temperature gluon medium, by adding an in- 
teraction of gluons witli a background Polyakov loop. It 
is well known that the latter is an exact order parameter 
for the confinement-deconfinement transition in the pure 
gauge theory (30.-i33ij ■ Therefore, having access to the 
Polyakov loop immediately leads us to the possibility to 
study the effect of the phase transition on the dynamical 
properties of gluons, mainly on their mass. The stan- 
dard quasiparticle approach, instead, tries to account for 
all the dynamics o nly by mean of temperature dependent 
masses 13"16|,[l^,[2l|,[2J|. This always leads to the neces- 
sity of diverging or steadily increasing masses as T —>■ Tc 
not only of the case of pure gauge theory, but also for the 
QCD case which includes quarks. We show that combin- 
ing a T-dependent quasiparticle mass, mg(T), with the 
Polyakov loop dynamics results in a quite different be- 
havior of TOg(T) as T ^ Tc. 

The combined approach of the quasiparticle picture 
with the Polyakov loop thermodynamics has been poorly 
investigated. To our knowledge, there is only a first pre- 
liminary study about the role of the Polyakov loop in 
the context of the quasiparticle description of the gluon 
plasma |l7l |. which has some overlap with our study. In 
the present article, we take a simple matrix model action 
as the effective potential for the Polyakov loop. Then, 
we use the Weiss mean field procedure to compute av- 
erages. Such a procedure has been used in the study of 
the QCD phase diagram by means of the Polyakov ex- 
tended Nambu- Jona Lasinio model 'sS, 'ss'] , see [s^, [s^l . 
In the next Section we describe the technical details, it 
suffices to say here that the main difference between this 
procedure and the one commonly used in the literature is 
that, in the former case, averages are computed in terms 
of integrals over the gauge group with the appropriate in- 
variant gauge measure having the pure gauge mean field 
action as the integration weight. We can anticipate that 
in such a way we are not only able to reproduce the large 
trace anomaly, Tj^ — e — 3P, with a quite good accu- 
racy, but also identify the different roles of quasiparticles 
masses and Polyakov loop in Yang-Mills thermodynam- 
ics. 

The paper is organized as follows. In Section II, we 
describe the matrix model in the Weiss mean field po- 
tential and its coupling to quasiparticle masses, showing 
the explicit derivation for SU{3) gauge group. In Sec- 
tion III, we compare the results of our approach to the 
IQCD thermodynamics discussing in particular the spe- 
cific contributions to the pressure and to the interaction 
measure of the Polyakov mean filed and of the quasipar- 
ticle masses . In Section IV we draw first conclusions of 
our work along with future possible investigations and 
developments of the present approach. 



Polyakov loop. The free energy of the model is expressed 
as a linear combination of two contributions: the first 
one describes the thermodynamics of the Polyakov loop; 
the second one, on the other hand, is the contribution of 
gluon quasiparticles coupled to the Polyakov loop. We 
specify the two terms of the effective potential below. 



A. The matrix model in the Weiss mean field 
potential 

In this Section we follow closely the notations of [H, 
[37j . The starting point is the action of a matrix 
model [Hi, 

Spm[L] = -N^C'/^Y^epixy^x + y) , (1) 



where x, y denote the lattice site and its nearest neigh- 
bors respectively. Moreover, 



(2) 



corresponds to the traced Polyakov loop in the funda- 
mental representation. The Polyakov line in the repre- 
sentation TZ of the gauge group is given by 



L-iiix) = T'exp 



I [ ' dT Al[x,T)Ta,n 
Jo 



(3) 



where Ta^-jz are the gauge group generators in the repre- 
sentation TZ. 

The model specified by Eq. ([Ij corresponds to a matrix 
model, the untraced Polyakov loops L(x) corresponding 
to the dynamical degrees of freedom. The partition func- 
tion can be written as 

Z = eM~^PM/T) = f'[ldL{x)eM-SpM[L]) . 

(4) 

Here, dL is the group invariant measure on the SU{Nc) 
gauge group. It can be written as 



dL= Y[ de^Y[\e" 

i—1 i<j 



-J, (5) 



and the integration over the variables 9i is performed over 
the range [0, 27r]. 

Adopting the Weiss mean field procedure, the action 
in Eq. ([T]) can be written as 



(6) 



II. THE MODEL 

The system we consider in this article consists of a gas 
of gluon quasiparticles, propagating in a background of 



where a, /3 denote the mean fields of the Polyakov loop 
and of its conjugate respectively. In terms of Smf, the 
average of an operator O, namely (O), is defined as 



{0) = 



fllxdL{x)0[L]c 
IUxdL{x)e-^ 



(7) 
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The thermodynamic potential f2 pm c&tl be written by 
adding and subtracting the mean field action Smf in the 



exponent of Eq. Q . We have 



log(e"(^^""^™^')) - log j J\dL{x)t 



where the average has to be performed as in Eq. ([7]). 



(8) 



Equation ([5]) is exact. In the mean field approximation, we replace (^e~^^^'^''~^'"f'') with e"^'^^""'^'"-'^ In this 
way, the thermodynamic potential depends exphcitely on the mean fields a, /3, and their values are determined 
selfconsistently by the stationarity condition. In the mean field approximation the thermodynamic potential can thus 
be written as [SGL \3T\ 



TV 



log dL e 



(9) 

In the previous equation, d corresponds to the spacetime dimension, so that 2{d — 1) is the coordination number of 
the cubic lattice, that is the number of nearest neighbors of the site x. Furthermore, is the lattice spacing, and V 
is the volume of the lattice; as a consequence, V/al is the number of lattice sites. 
The thermal average of the Polyakov loop and of its conjugate are given by 



JdL 



[cos(A^c/33^f)3?^f - sin(iVc/33^F)3^/ 



/ dL e^-"^^" cos{NcP'^£f) 
J dL e^""^^^ [cos(7Ve/3S£F)3^^^F + sin(jVe/3S£f 
/ dL e^"'^^''^'^ cos(7Vc/?3^f) 



(10) 
(11) 



From the above equations it is then clear that {£p) and {£*p) are different when /3 ^ 0. This happens for example in 
theories with fermions at finite baryon chemical potential. 

In the following we focus on the case of the pure gauge theory, in which {£p) = {£*p). We have thus (3 — 0, and 
Eq. ^ can be written as 



PM 



bT 



{-2{d - l)Nle-"''^{£pf + [a{£F)] 



log I dL e 



NJamF 



(12) 



where we have used the shorthand notation b = V/a^. In the following, a and b in the above equation are treated as free 
parameters, which will be fixed, together with the remaining two parameters of the model which will be described in 
the next Section, by a minimization procedure of the mean quadratic deviation between the theoretical computations 
and the Lattice data or pressure, energy density and interaction measure. For the sake of future reference, we call the 
model described by the above thermodynamic potential as the pure matrix model (PMM). 

I 



B. The quasiparticle gas contribution 

The gluon quasiparticle contribution to the thermody- 
namic potential reads [l^ 



2T 



d^k 
(2^ 



Tr^ log 1 - La e 



-E(k)/T 



(13) 



In the above equation. La corresponds to the Polyakov 
line in the adjoint representation, as defined in Eq. (O, 
and the trace is taken over the indices of the adjoint rep- 
resentation of the gauge group. M oreover, t he quasipar- 
ticle energy is given by E{k) = y/W + AP, where M is 
supposed to arise from non-perturbative medium effects. 



In the case of very high temperature, where the gauge 
theory is in the perturbative regime and thus M = 0, it 
can be proved that Eq. (IT5t is the effective potential for 
the adjoint Polyakov loop [33|- At lower temperatures, 
where non-perturbative effects are important, Eq. (jl3p 
must be postulated as a starting point for a phenomeno- 
logical description of the thermodynamics of the gluon 
plasma. In this article, we make use of a temperature 
depentent mass whose analytic form is given by (20l | 



Stt^T^ 



11 log [X{T~w)r 



(14) 



The parameters w,X are determined by requiring that, 
for pressure and energy density, the mean quadratic devi- 
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ation between our theoretical computation and the Lat- 
tice data of is a minimum. In particular, we find 
w = 240 MeV and A = 0.25, see Section III. The tempera- 
ture dependent quasi-particle mass in Eq. (|14p has several 
motivations: firstly, in the very large temperature regime 
in which QCD is pcrturbative, one has AI — g{T)T; 
moreover, for T close to Tc there are remnants of the 
confining mechanism and, a priori, there is no reason 
to reduce these effects to the Polyakov loop dynamics; 
finally, for T — )• T~, the glueball mass decreases [25l |. 
Mgiuebaii ~ 1 GeV at T = Tc , and in a valence gluon 
model one expects an effective mass of the quasiparticle 
M « 500 MeV at the melting temperature. 

In order to have a combined description of the gluon 
plasma around the critical temperature, in terms of the 
Polyakov loop (which acts as a background field) and of 
gluons quasiparticles (which propagate in the Polyakov 
loop background) , we add Eq. (jl3l) to the strong coupling 
inspired potential in Eq. (|12p . In order to do this, and 
to be consistent at the same time with the Weiss mean 
field procedure outlined above, we follow the lines of 
and write the thermodynamic potential as 



PM 



(15) 



where the average must be understood as defined in 
Eq. (O. In taking the average, we use the same approx- 
imation of (36| and write 

logdet (l - La e^^/^) ^ = log (|det (l - La e-^/^ 

(16) 

This approximation is also discussed in the Appendix 
B of [331, where it is proved that taking the average of 
the log (that is, {^qp)) instead of the argument of the 
log, would lead to unphysical artifacts in the case the 
particles are taken in representations different from the 
fundamental one. We will discuss this aspect in more 
detail in the next Section. 



C. The SU{3) case 

In the case of the color group SU{3) the Polyakov 
line can be parametrized in terms of two angles, 63 and 



6*8, namely Lp = diag(e**'^ e'^'s, e 



-i{e3+98) 



); the invariant 



J 



As a result we find 

log (^detA ( 1 
where we have defined 



Haar measure reads: 



dL = 



67r2 



[sin(6'3 - 6*8) - sin( 



3) + sin( 



(17) 



with Oq = -{63 + 9s). We define 

Fia) = JdL e"^^^^- = J dLe^'^f^'^''^'^ (18) 

where 

Moreover we have 

1 /dLe"*^/(»i'''^)/(^i,02) 



(19) 



Nr. 



F{a) 



(20) 



The thermodynamic potential for the pure matrix model 
in this case then reads 

^PM = hT {-QNle-"''^ {tpf + Nr aiip) - logF(a)) . 

(21) 

The quasiparticle contribution for the SU (3) case in 
the mean field approximation reads 



^qp 2T 



d^k 
(2^ 



log ( dct A ll- La e 



-E/T 



, (22) 



where the determinant is taken in the adjoint color rep- 
resentation, and the adjoint Polyakov line is given by 
La = exp [i {03Tl + esTl)], with (r^)ac the appropri- 
ate generators in the adjoint representation. We make 
use of the following equation to express the adjoint loop 
in terms of the fundamental one: 



L'X = 2Tr 



LpTaLpTi, 



(23) 



where Ta, Ti, correspond to the generators of SU{3) in the 
fundamental representation, normalized with the condi- 
tion TT[TaTb\ — Sab/'^- To compute the functional deter- 
minant in Eq. (1221) . we chose a representation in which 
La is diagonal; the eigenvalues of La are Ai = A2 = 1, 
A3 = A^ = e'(^^-^»\ As = A* = e*(''^-''«), and finally 

A7 = a; = e*(««-s9), with Og = -{e3 + es). 



La e 



-E/T 



21og(l 



3 

E 



log 1 



2 (w.) e 



-E/T 



Wi = COS (6*3 - 6*8) , 
W2 = cos (6*3 26*8) , 
a;3 = COS (26*3 + 6*8) , 



(24) 



(25) 
(26) 
(27) 



and the (...) denotes the average over the gauge group elements as defined as in Eq. ([T]). 



5 



Our expression for the functional determinant is in 
agreement with Eq. (4.5) of [s?!, once we put x = e~^^'^ 
and change x —x in Eq. (j22l) to take into account that 
in [37j . fermions in the adjoint representation are con- 
sidered instead of bosons. It can be proved anahtically 
that 



1 



a = 



moreover, numericahy we have checked that 



^1 



^2 



1 



q: — )■ oo 



(28) 



(29) 



For the sake of reference, we write down the equation 
for the traced fundamental and adjoint Polyakov loops in 
terms of the angles 6*3 and 6*8, namely 



1 + LOl + UJ2 + LU3 



(30) 



63 — 6s O3 + 29s 26*3 + 8s , 

I A = COS COS COS . (31) 



We have checked numerically that {If), {^a) ~ 1 in the 
asymptotic limit a — >■ 00 and {£f), {^a) — for a — > 
(the latter step can be proved analitically as well). The 
asymptotic behavior of the averages of the uJi in Eq. 
implies that in the deconfinement phase, at very large 
temperature where {£a) ~ 1, the thermal distribution 
of the quasigluons approaches that of a perfect gas of 
massive particles. However the thermodynamics of the 
system remains different from the one of a mere massive 
gas because of the Polyakov background mean field. 



III. NUMERICAL RESULTS 

In Fig. [T]we plot the pressure p = —il (red data) and 
energy density (green data), normalized to the Stefan- 
Boltzmann value, as a function of temperature (mea- 
sured in units of the critical temperature), for the SU(3) 
case. Given the pressure, the energy density is com- 
puted by virtue of the thermodynamical relation per- 
forming the pertinent total derivative of the tempera- 
ture, e = —p + Tdp/dT. Of the four parameters in our 
model, one of them is fixed in order to reproduce the 
first order phase transition at T = 270 MeV; the re- 
maining three parameters are fixed in order to require 
that the mean quadratic deviation for pressure, energy 
density and interaction measure between our computa- 
tion and the Lattice data is minimized. This procedure 
leads to the numerical values a = 901.9 MeV, b — (157.5 
MeV)3, w = 240 MeV and finally A = 0.25. These pa- 
rameters produce a gluon mass M — 360 MeV at T = Tc. 
In the figure, psB — (iA7''^r^/45 and €sb — 3pss where 
dA = -I is the dimension of the adjoint representa- 
tion. Dots correspond to lattice data taken from Ref. Q ; 
solid lines are the result of our numerical computation 
within the quasiparticle model. 




■ energy density (model) 

■ energy density (lattice) 
- pressure (model) 

pressure (lattice) 
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FIG. 1. Upper panel. Pressure (indigo, solid line) and energy 
density (orange, dashed line), normalized to the respective 
Stefan-Boltzmann values, as a function of temperature (mea- 
sured in units of the critical temperature Tc = 270 MeV), for 
the SU(3) case. Dots correspond to lattice data taken from 
Ref. solid lines are the result of our numerical computa- 
tion within the quasiparticle model. Lower panel. Interaction 
measure per degree of freedon as a function of temperature. 
The dashed orange line corresponds to the result of Ref. 



In the lower panel of Fig. [T] we compare our results for 
the interaction measure, A — {e — 3p)/T*, with the lat- 
tice data (represented by dots in the figure). We notice 
that our model reproduces fairly well the peak of the in- 
teraction measure observed on the lattice in the critical 
region. We also compare our results with the ones ob- 
tained in [TtI , where a phenomenological potential for the 
Polyakov loop is considered instead of our matrix model, 
and the usual mean field approximation is used. The 
comparison shows that the peak is much better repro- 
duced in our case. However one could argue that the tail 
at high T in instead not so well reproduced; in the next 
section we will see that such a drawback can be cured by 
a slightly different Mg{T) behavior. 

In the upper panel of Fig. [5] we plot our result for 
the gluon mass a gain st temperature, and we compare it 
with the result of I18| obtained without the Polyakov loop 
dynamics. In Ref. [l8| the following ansatz is used: 



M = 



Bt , 



(32) 



{t-5Y 

with t — T /Tc. Numerical values of the parameters are 
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FIG. 2. Upper panel. Ratio AI/T for our model (indigo, solid 
line) and for the model without Polyakov loop (orange, dashed 
line) of 18], against temperature, the latter measured in units 
of the critical temperature T = 270 MeV. Lower panel. Ex- 
pectation value of the fundamental (indigo, solid line) and 
adjoint (orange, dashed line) loops against temperature. 



A = 1.42 Tc, B = 0.533 Tc, c = 0.46 and S = 0.952. The 
first addendum on the r.h.s. of the above equation, which 
is interpreted as the energy contained in the correlation 
volume, is important close to Tc, and it is necessary to 
have a large mass as T approaches Tc from above. In 
fact, when the Polyakov loop is not introduced, a large 
mass is needed in order to suppress pressure and energy 
density around Tc. In the case the phase transition is 
of the second order, as it happens for SU{2), correlation 
length diverges at Tc] therefore, (5 = 1 since the mass 
diverges as well. When the phase transition is of the 
first order, the correlation lenght (hence the correlation 
volume) does not diverge at Tc, but it is expected to be 
large; as a consequence, the mass should increase as T 
approaches Tc from above. For this reason S in Eq. ([5^ 
is slightly less than unity ^ < 1 . The second addendum 
on the r.h.s. of Eq. dominates the mass at large 
temperature, and its origin is of perturbative nature. 

In our case, the statistical suppression of states be- 
low and around Tc is achieved by virtue of the Polyakov 
loop. For this reason, we do not need to have a large 
mass as T approaches Tc. This mechanism is similar to 
the statistical confinement that takes place in the PNJL 
model [13, |4l| . To show this aspect of the model clearly, 
we have plotted in Fig.[3]the full pressure obtained within 




full pressure 
qp and Polyalovloop 
pure qp gas 



1.0 1.5 2.0 2.5 3.0 3.5 



full result 
pure qp gas 
without dm/dT 




FIG. 3. Upper panel. Full pressure of our model (indigo, 
solid line), perfect gas pressure (orange, dashed line) and 
quasiparticle contribution (teal, dot dashed line) against tem- 
perature. All the data are measured in units of the Stefan- 
Boltzmann pressure. Lower panel. Interaction measure of our 
model (solid, indigo line) compared with the pure quasiparti- 
cle model (orange, dashed line). In the figure, the dot dashed 
teal line corresponds to the interaction measure we would ob- 
tain neglecting the derivative of the gluon mass, see the text 
for more detail. 



our model (indigo solid line), and we have compared this 
result with the pressure of a perfect gas having the tem- 
perature dependent mass we have obtained in our com- 
putations (orange dashed line). From the data shown in 
the figure, the role of the Polyakov loop to suppress ther- 
modynamically the contribution of states is noticeable. 
For completeness, we have also plotted the quasiparticle 
gas contribution to the pressure, ^qp, with Polyakov loop 
taken into account (the latter has been selfconsistently 
computed at any temperature). 

The role of the Polyakov loop is apparent also in the 
lower panel of Fig. [3l where we compare our results for 
the interaction measure, which has Polyakov loop ther- 
modynamics built into it, and the interaction measure for 
a gas of quasiparticle with mass given by Eq. (jl4D . For a 
matter of simplicity, the latter is estimated by using the 
relativistic Boltzmann distribution function for gluons, 
namely [2^ 



Ab = 2 



(7V2 _ 1)M2 



'M 


/dMV 







Ki{M/T) , (33) 



27r2r2 

the analytical expression allows to identify the two con- 
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T/r, 

FIG. 4. Upper panel. Normalized pressure and energy density 
against temperature. Middle panel. Interaction measure per 
degree of freedon as a function of temperature. Lower panel. 
The ratio M/T against temperature. 



tributions to As, we have verified that turning to the 
Bose-Einstein statistics does not lead to significant dif- 
ferent results. 

It is evident that the Polyakov loop gives a relevant 
contribution to the interaction measure. Moreover, it is 
instructive to show the concrete effect of a temperature 
dependent mass on the interaction measure. To this end, 
in the figure we have plotted Eq. in which we neglect 
the term proportional to dM/dT: the result is shown by 
the teal line. Since dM/dT in our case is positive, its 
contribution suppresses the interaction measure, as it is 
clear from the data shown in the figure. 

As a final investigation, we have considered another 



ansatz for the thermal gluon mass, that is 

Munear - Mo + h (^^ - 1^ , (34) 

where Mq = 500 MeV and h = 126 MeV (these two co- 
efficients have been fixed once again by a minimization 
procedure of the quadratic relative deviations). The re- 
sults are summarized in Fig. |4l Few comments follow. 
Secondly, the peak of the interaction measure is slightly 
suppressed: in this case the suppression is of about 10%, 
which has to be compared with the case of the logharith- 
mic mass (less than 2%). However the tail of A is bet- 
ter reproduced, because of the combined effect of mass, 
Polyakov loop and dM/dT. The role of the Polyakov 
loop in avoiding the divergence of the gluon mass as Tc 
is approached from above remains. We leave a more sys- 
tematic study of the several possible ansatzes to a future 
work. 



IV. CONCLUSIONS 

In this article, we have studied a combined description 
of the deconfinement phase of the SU{3) gluon plasma, in 
terms of gluon quasiparticles and Polyakov loop. The two 
are mixed by a phenomenological potential, see Eq. (fTB]) , 
which is inspired by weak coupling computations [39] . To 
the quasiparticle contribution, we have added a potential 
for the Polyakov loop given by Eq. (TT^ ; its analytic form 
is inspired by the leading order strong coupling expan- 
sion, and it has been used in the literature in the context 
of the PNJL model 

Our main purpose is twofold. Firstly, we are inter- 
ested to a simple description of the lattice data about 
the thermodynamics of the gluon plasma. This is in- 
teresting because it allows to understand which are the 
relevant degrees of freedom in the deconfinement phase 
of the theory. Moreover, once the proper degrees of free- 
dom are identified, the effective description studied here 
can be completed by adding dynamical quarks. Employ- 
ing a relativistic transport theory, this can allow also a 
direct connection between the developments of effective 
models and the study of the quark-gluon plasma created 
in relativistic heavy-ion collisions. 

Within our model, we are able to reproduce fairly well 
the lattice data about the gluon plasma thermodynamics 
in the critical region. In particular, we reproduce the 
peak of the interaction measure within a 2%, see lower 
panel in Fig. [1] 

One of the main conclusions of our work is that because 
of the coupling to the Polyakov loop, a gluon mass of the 
order of the critical temperature, rrig ~ 1 — 2 Tc, is enough 
to reproduce the lattice data in the critical region. This 
is different from what is found in standard quasi-particle 
models. The Polyakov loop is neglected and the gluon 
plasma is described in terms of a perfect gas of massive 
gluons. As a consequence, a mass rapidly increasing with 
T is needed to suppress pressure and energy density in 
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the critical region. On the other hand, in the case studied 
in this article, the suppression of the pressure and energy 
density is mainly caused by the Polyakov loop, as shown 
in Fig. [21 Therefore, lighter quasiparticles with a smooth 
T— dependence describe fairly well the lattice data in the 
critical region. 

There are several interesting directions to follow to ex- 
tend the present work. Firstly, having in mind the study 
of the quark-gluon plasma, we plan to add dynamical 
massive quarks to the picture. It is also of a certain in- 



terest to span the parameter space in more detail, eventu- 
ally analyzing several functional forms of the gluon mass. 
Moreover, even in the case of the pure gauge theory, it 
would be interesting to investigate the way to adjust the 
interaction measure in order to reproduce lattice data 
in the regime of very large temperature. Even more, it 
would be of a certain interest to extend our study to the 
case of different gauge groups. 

Note added. While preparing the present article, 
Ref. [13 was submitted to arXiv, which has some overlap 
with our work. 
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